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ABSTRACT 


The work presented is a preliminary feasibility study of strato- 
spheric and mesospheric water vapor measurements using the micro- 
wave lines at 22 GHz (22.235 GHz) and 183 GHz (183.31 GHz). 

The resonant cross sections for both the 22 GHz and the 183 GHz 
lines are presented and used to model the optical depth of atmospheric 
water vapor. The range of optical depths seen by a microwave radio- 
meter through the earth's limb is determined from radiative transfer 
theory. Radiometer sensitivity, derived from signal theory, is 
compared with calculated optical depths to determine the maximum 
height to which water vapor can be measured using the following methods: 
passive emission; passive absorption; and active absorption. 

The report concludes that measurements using the 22 GHz line 
are limited to about 50 km whereas the 183 GHz line will enable 
measurements up to and above 100 km for water vapor mixing ratios 
as low as 0. 1 ppm under optimum conditions. 



CHAPTER I 


INTRODUCTION 


The question of how much water vapor exists in the stratosphere 
and above is one that has received attention for a considerable number 
of years but which has not yet bee” fully answered. The study of water 
vapor concentrations is important since water vapor and its dissociated 
products are active in the photochemistry of the stratosphere and meso- 
sphere and thus effects ozone and contributes to the escape of hydrogen, 
for example. In addition it is thought to be responsible for hydrated 
ions, OH air glow emissions, and may control ionization loss processes 
as high as 80 km. Thus the amount of water vapor in the upper atmo- 
sphere is an important problem since it plays an important part in the 
chemistry and heat budget of the atmosphere. 

On the basis of the above, it would appear that the measurement 
of water vapor concentrations in the upper atmosphere is of considerable 
importance to the scientific community. The best measurements of 
water vapor below 30 km were made with frost point hygrometers; how- 
ever, these do not appear to be sensitive enough to the small concentra- 
tions of water vapor above about 3 0 km. They also suffer from the 
problem of local water vapor contamination. 

The possibility of using radiometric methods for the detection 
and measurement of water vapor in the stratosphere and mesosphere 
was considered in a recent report (Longbothum, 1974). That report 
discusses the probable concentrations of water vapor and what is 
known about the microwave, millimeter wave, and infrared regions of 
the water vapor spectrum. The report concludes with a discussion of 
various methods of radiometric detection. The next logical step in 


evaluation of a radiometric detection system is to examine the Edinimran 
sensitivity of the various systems. This is the purpose of the present 


CHAPTER. II 


ABSORPTION CROSS SECTIONS 


The microwave and millimeter wave regions will be examined 
in this report (Moody, (1971) reported that these were the only bands 
suitable for making measurements of water vapor concentration above 
30 km). There are only two resonant water vapor lines in the micro- 
wave and millimeter wave regions: one at 22.235 GHz (1.35 cm); and 
the other located at 183.31 GHz (1.64 mm). The resonant cross 
sections are computed from relations A. 1, A. 2, A. 3, A. 4, and A. 6 
in appendix A using the relation 


o- (h) 


Ka (h) 
N (h) 


( 1 ) 


Table 1 is obtained for the centers of the resonant lines at 22.235 GHz 
amd 183.31 GHz. The above cross sections are plotted in Figure 1. 

The plot shows that the cross sections can be described (for the altitude 
range 3 0-80 km) by the function 


tr (z) = cr e 
o 


z/C 


where 


( 2 ) 


_ Az 

A(lntr(z)) 


C 


7.4 km 


HEIGHT IN KM 



Cross Sections 



Table 1; Absorption Gross Sections 


Av 


22 


Av 


183 


°"22 (h) 


^183 ^ 


30 km 

4. 0 

X 

I0 7 

Hz 

4. 0 

X 

io 7 

Hz 

1. 0 

X 

10- 22 

2 

cm 

2.6 

X 

io - 20 

2 

cm 

40 km 

9.7 

X 

10 6 

Hz 

9.7 

X 

io 6 

Hz 

4.0 

X 

io ' 22 

9 

cm" 

9.5 

X 

io - 20 

2 

cm 

50 km 

2.7 

X 

10 6 

Hz 

2.7 

X 

io 6 

Hz 

1.5 

X 

10" 21 

2 

cm 

3.1 

X 

H>- 18 

2 

cm 

60 Ion 

6« 6 

X 

10 5 

Hz 

7.1 

X 

io 5 

Hz 

6. 1 

X 

io - 21 

2 

cm 

1.3 

X 

io ' 18 

2 

cm 

70 1cm 

2,3 

X 

10 5 

Hz 

3.3 

X 

10 5 

Hz 

1.7 

X 

10- 20 

2 

cm 

3.9 

X 

!0- 18 

2 

cm 

80 km 

5. 0 

X 

4 

10* 

Hz 

2„1 

X 

10 5 

Hz 

6.6 

X 

10“ 20 

2 

cm 

8.3 

X 

io - 18 

2 

cm 

90 Ian 

2. 6 

X 

10 4 

Hz 

2. 1 

X 

10 5 

Hz 

1.3 

X 

10“ 19 

2 

cm 

7.8 

X 

IO ' 18 

2 

cm 

100 1cm 

2.8 

X 

10 4 

Hz 

2.2 

X 

10 5 

Hz 

1.4 

X 

io " 19 

2 

cm 

5.8 

X 

IO ' 18 

2 

cm 

110 1cm 

3.0 

X 

10 4 

Hz 

2-5 
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io 5 

Hz 

1.3 

X 

io - 19 

2 

cm 

3.7 

X 

ID" 18 

2 

cm 

120 km 

3.5 
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10 4 

Hz 

2.9 

X 

10 5 

Hz 

1.0 

X 

IO " 19 

2 

cm 

1.9 

X 

IO ' 18 

2 

cm 


* (See Appendix A) 



CHAPTER HI 


OPTICAL DEPTH OF WATER VAPOR ABOVE 30 KM 


The transmission of a plane atmosphere is described by the 
optical depth t (v, h, (p) defined as 


t (v, h, <p) - ^ tr ^ (z, v )n^ (z) sec <p dz (3) 

i h 


where 
cr 1 (z, v ) = 


o- 


scatter 


(z, V ) + cr- 


abs orptlon 


(z, v ) 


2 

- the total extinction cross section, in cm , for the ith 
constituent 

i -3 

n (z) = the number density of the ith constituent in molecules ° cm 

<j> = the zenith angle in degrees 

Spherical atmospheres can be treated with this expression for 
zenith angles less than about 80°, However, for zenith angles greater 
than 80°, the Chapman Function must be used in place of sec cp 
{Rishbeth, et al. , 1969)= For the sake of simplicity, the plane atmo- 
sphere solution will be used for this report- 

For the molecules, atoms, and particles above 30 km, 
non- resonant scattering can be neglected in general. Thus only the • 
resonant absorption cross section will be retained in the expression 
for optical depth. Also in the frequency spectrum near the two water 
vapor lines, only four other molecules have resonant absorption 
(emission) lines (Longbothum, 1974). They are O^, SO^, NO^, and 
N^O. The concentrations of SO^, NO^ and N^O are thought to be 
very small. The concentration of above 30 km, follows quite 
closely that of H^O and must be considered in the case of the 183.31 GHz 
water vapor line. This possible complication will be dealt with in a 



7 


future study. For the present only water vapor will be considered. 

It was suggested in the recent report (Longbothum, 1974) that 

the most probable water vapor concentrations lie between 0. 1 ppm and 

10,0 ppm. Thus a value of n corresponding to 0. 1 ppm and a value 

of n corresponding to 10 ppm will be used to calculate T . and 
max mm 

r respectively^ 
max * 3 

The product cr(z)n(z), the differential optical depth, is plotted 
in Figure 2. This shows that for the height range up to about 80 km, 
the product tr (z) n (z) is relatively constant and approximately equal to 
o-^n^. Above 80 km, the product is height dependent and can be 
described by the function 

<r(z)n(z) = o'qELq exp [= (z - 80 km) / H 1 ] (4) 


tt) ** As __ 

Ain l_n(z)° <r (z)j 
for the 22.235 GHz line 


and 


H» = 5 km for the 183.31 GHz line. 

The vertical optical depth ( <p = 0) can be integrated for the two 
ranges of z to yield 

f n Q cr o [ (80 km - h) + H'] h < 80 km (5) 

| & o <r Q H* exp [-(h = 80 km)/H'] h > 80 km 

_iq „ i 

The average value of n Q cr is 3.4 x 10 cm for the 22.235 GHz 

— O a 1 

line and 8. 0 x 10 cm" for the 183.31 GHz line. These values are 
based on Figure 2. Table 2, also obtained using these values, shows 
that attenfcuation will be slight at 22.235 GHz and quite significant at 
183.31 GHz. 



HEtGHT IN KM 



30 km 1. 

40 km 1. 

50 km 1. 

60 km 8. 

70 Ion 5. 

80 km 2. 

90 km 3. 
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110 km 1. 

120 km 2. 
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54 x 10“ 5 

1.54 x 10" 3 

20 x 10" 5 

1.20 x 10“ 3 

70 x 10" 6 

8.70 x 10" 4 

35 x lO" 6 

5.35 x 10" 4 

01 x 10" 6 

2.01 x 10" 4 

79 x 10" 7 

3.79 xlO" 5 

16 x 10" 8 

7. 16 x 10“^ 

35 x 10" 8 

1.35 x 10" 6 

56 x 10" 9 

2.56 x 10“ 7 


4.38 x 10 


4.38 x 10 


3.59 x 10" 3 

3.59 x 10’ 

2.79 x lO' 3 

2. 79 x 10 

1.99 x lO" 3 

1,99 x 10 

1.20 x 10" 3 

1. 20 x 10 

3.98 x 10" 4 

3.98 x 10' 

5.39 x 10“ S 

5.39 x 10' 

7.30 x 10" 6 

7.30 x 10' 

9.88 x 10" 7 

9.88 x 10' 

1 . 34 x 10" 7 

1 . 34 x 10 



CHAPTER XV 


RADIATIVE TRANSFER 

To measure water vapor concentrations using a radiometer, one 
must establish the relationship between various sources of radiation, 
the atmospheric optical depth t, and the radiation that the receiving 
antenna intercepts. This relation is given by the theory of radiative 
transfer. 

Consider an antenna connected to a receiver. The antenna is 
viewing a source (temperature T^) at a great distance through the 
relatively cold atmosphere (temperature Tg) that lies between the two. 

This atmosphere has an optical depth t(v o ). If one includes the atmo- 
sphere which lies beyond the source also, then the total optical depth 
is t' ( v^). The temperature that the antenna will see then is given by 
the solution of the equation of radiative transfer (Barrett, 1964). 

T ( V Q ) = T A e " T( ‘' V o ) + T b (6) 

where 

T (v ) is the received antenna temperature for frequency v q 

T A is the unattenuated -temperature for any source, finite or 
extended in the antenna beam. 

Tg is the excitation temperature of the gas (background temperature), 

t ( v ) is the optical depth of the gas between the observer and the 
source 

and 

t' ( v q ) is the optical depth of the total gas including any beyond the 
source. 

When the radiometer is tuned to v (not a resonant frequency of 
the atmospheric gas), the optical depths change to new values t(v ) and 
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T* (v) and the change in antenna temperature between T (v q ) and T (v) 
is 


AT = T(v q ) - T (v) = T a [e 


- T < v o l . e ' Ttv) ] 


+ T b [e 


-T' (V) -T' ( V ) 


- e 


V] 


( 7 ) 


AT represents the difference signal measured by a receiver as it is 
tuned through the resonant line atv^, However, AT will include both 
absorption and emission effects of the gas. In order to make a more 
accurate measurement, there must be some method of separating the 
two responses. If the antenna is directed away from the source or 
(in case of an active experiment) the transmitter is turned off, then A T 
is the signal from the emissions of the gas. 


AT 


emission 


= AT' (v Q> v) 



[e 


- T» ( V ) 



£ T-o T' (V ) for 1 > > T' (V ) > > T* (v ) 
B o o 


( 8 ) 


This approximation is generally valid since the values r' (i> ) and r(v) 
are due to all the nonresonant extinction processes such as wings of 
higher order lines and scattering. Thus they are much smaller than 
the resonant optical depths r 1 (v^) and t(v q ). 

If AT' is subtracted from A T, the contribution due to absorption 
only is obtained. 


AT 


absorption 


AT" (v ,v) = T. [e“ T( V - e" T(v) ] 

O A 


H -T. t(v ) for 1 > > t(v ) > > t(v) 
A o O 


( 9 ) 
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The measurement only involves the magnitude of the signal thus the 
relations between the receiver input and the atmosphere are 


I AT' (v Q ) |S T b t* (v o ) 
for an emission experiment and 


( 10 ) 


| AT" <V Q ) | £ T a t (v q ) (11) 

for an absorption experiment. 

For most conceivable experiments 

t (v ) s T 1 ( v ) 

' o o 

so that a common expression can be written for both measurements 
| AT{* o )|£ Tg r { v o ) (12) 

where Tg is the source temperature and is given by 


r T^ for absorption 

Tg — j 

Tg for emission 

This is the equation that couples the atmosphere to the radiometer. 
However, in order to detect a water vapor concentration, one must 
have 


I AT <v o ) 


> AT . 


min 


or equivalently 


I AP (v Q ) 



(V o } 


> AP 


min 


(13) 
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where AT . is the minimum detectable temperature the radiometer 
mm 

can sense and is limited by the noise in both source and receiver; 

AP . is the minimum dete ctable change in signal power that can be 
mm 

sensed by the radiometer and Pg is the radiated power of the source. 


CHAPTER V 


T SN " T R + T L 


RADIOMETER SENSITIVITY 
The expression for AT^ (see appendixes B and C) is 

AT . = 2T SN” 

mm — 

\/ B fe (B. 24) 

for a voltage signal to noise ratio a and where 

Tgj^ = T^ + Tg + T^ for a thermal source and 

T pm - T_ + T t for a nonthermal source. 

SN R L> 

T^ = the noise temperature of the receiver 

T_ = the noise from antenna losses, sidelobe noise, etc. 

J-J 

B = the receiver bandwidth in H 

z 

and t is the post detection integration time in seconds of the receiver. 
Notice that if Tg > T^ + T^ then there is nothing to be gained by- 
building a better receiver since the noise fluctuations will be due to the 
source. Thus the limiting factor on absorption measurements will be 
the source noise in the case of a thermal source. 

A minimum measurement can be made only if the equality holds 
in equation (13). This measurement can be enhanced by slant path 
viewing, thus the relationship becomes 


t cy 

sec <p 


mm 
sec <p 


AT . 
mm 

Tg sec 0 


2 

nTbT 


a 

sec <p 


t r + 



15 


It would appear that one could reduce T by increasing B or t in 

the above equation. However, t is limited by the amount of available 

observing time and desired spatial resolution. This is determined by a 

change in source position as a function of time, by a time dependent 

change in the power received by the antenna, or by the time of flight of 

the vehicle. Also, t can. not be increased beyond the point where it 

begins to distort the true source profile (Kraus, 1966). For rockets, the 

0 T 

time of flight would limit the integration time to 10 - 10 seconds. For 

1 2 

a satellite, the time range would be about 10 - 10 seconds. For a 

1 3 

balloon, the time range would be about 10 - 10 seconds. 

The bandwidth B can not be made too wide or there will be a loss 
of spectral information and interference from radio signals of terres- 
trial origin (Kraus, 1966). Also B will be limited by the bandwidth of 
the signal, the current state of equipment development, and many other 
factors. B has a roaximum value of about 10^ Hz, but is only about 
200 MHz for low noise receivers (Staelin, 1969). 

System noise temperatures of radiometers have for years been 
the limiting factor in sensitivity. Values of system temperature have 
been about 3000°K for 22. 235 GHz radiometers to about 70, 000°K for 
183.31 GHz radiometers, (See Figure 3). This figure shows the state 
of the art system noise temperatures that existed 5 to 10 years ago. 
However, there has been a major breakthrough in mixer development in 
the last few years. Current values for system noise temperatures 
range from about 600°K for 1.3 cm radiometers to about 4000°K for 
2 mm radiometers (Buhl et al., 1971). Staelin (1973) also reports that 
system noise temperatures less than 300°K at frequencies up to 200 GHz 
may he possible with cooled Schottky barrier diodes. Waters (1974) 



NOISE TEMPERATURE IN °K 
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Figure 3 Noise Temperatures of Receivers 
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gives a state of the art system temperature of 1Z00°K at 183 GHz. Thus 

for this report the value of T R + will be taken as 600°K for the 

22. 235 GHz line and 1200°K for the 183.31 GHz line. 

The apparent source temperature, Tg, must be specified for 

each possible experiment. For atmospheric emission T~ is assumed 

to be 300°K, For an absorption experiment using the solar disc as an 

external source, T g S 10, 000°K for the 22.235 GHz region and 

Tg £ 6, 000°K for the 183.31 GHz region of the solar spectrum 

(Linsky, 1973). An active experiment using a 10 milliwatt broadband 

source (for a receiving antenna effective area A g = 10 M , a trans- 

ft 

mitter antenna gain of 100, and a distance of 3 x 10° meters between 
the receiver and transmitter) has a value of Pg 


p s = 


GA 


4 tr R 


8.84x10 ^ Watts 


(15) 


CHAPTER VI 


COMPARISONS OF EXPERIMENTAL CONFIGURATIONS 

There are three basic experiments that can be performed now. 
An emission experiment, a passive absorption experiment, and an 
active absorption experiment. These can be performed at 22.235 GHz 
or at 183.31 GHz. The experiment at 22,235 GHz suffers from weak 
line strength, but has the advantage of greater receiver sensitivity 
and less interference from other lines. The experiment at 183.31 GHz 
has a resonant cross section on the order of 200 times greater, but 
suffers from less sensitive detectors, and more interference from 
higher order lines so that equipment having very good resolution will 
be required. Both lines will be considered below, six figures are 
presented corresponding to the six possibilities to follow, B will be 
taken to be Av of the line and sec <£will be taken equal to 10. This 
will correspond to an occultation experiment (Longbothum, 1974). 

6.1 Passive Resonant Emission 

This experiment has the advantage that it does not require a 
source and can be performed at any time. The method is, however, 
an order of magnitude less sensitive to water vapor concentrations 
(Longbothum, 1974). 

The minimum optical depth that can be detected is given by 

AT min _ 2a I 1 60 \ 

T min - T g sec = tj rK - • I 10 T g 

( 16 ) 


for the 22.235 GHz line, and 
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min 


AT 


min 
sec $ 


2 a 

•itAv 




(17) 


for the 183,31 GHz line. 

The values of Tg in these expressions is taken as the average of the 
temperatures at height h and at 5 km above h. The values of T c 
obtained are (U. S. Standard Atmosphere, 1962). 


h 


h 

c 

30 km 

237 °K 

80 km 

181°K 

40 km 

2S7°K 

90 km 

1 88°K 

50 km 

269°K 

100 km 

224°K 

60 km 

248°K 

110 km 

286°K 

70 km 

210°K 

120 km 

41 1°K 


The maximum signal to noise ratio a r is taken as 10 and the minimum 
is taken as 1. The maximum time t is taken as 100 sec and the 
minimum time is taken as 1 sec. Thus the ratio ar/^yy ran S e from 

10 to 0. 1. 

Table 3 of minimum detectable optical depths is obtained from 
equations 16 and 17. These values reflect the increased ratiometer 
sensitivity at lower frequencies (see figure 3). 

Based on the limiting value of atmospheric optical depth 
presented in Table 2 and the minimum detectable optical depths 
presented in Table 3, Figures 4 and 5 are obtained. These figures 
show that at 22.235 GHz, water vapor can be measured to about 70 km 
under optimum conditions; However, the more realistic limit would 
be about 50 km. When the experiment is performed at 183,31 GHz, 


Table 3: Minimum Detectable Optical Depths for Passive Emission ( AT ^^ /Tg. sec $ 
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the situation is improved. For optimum conditions, measurements 
can be made up to 100 km with a more realistic limit being 80 km. It 
is clear from the graphs that an observation time of at least 10 
seconds should be used. 


6. 2 Passive Resonant Absorption 

This method is quite. sensitive but requires' a natural source. 

The strongest natural source is the solar disc. However, this method 
will be limited by the position of the sun in relation to a moving vehicle. 
This will impose short observational times on the measurement. The 
method does have the advantage that only one vehicle is required. 

The minimum detectable optical depth is 


T . 
min 


AT 


mm 


Tg sec (p 


2 a 

nT t A v 


for the 22.235 GHz line, and 


J- + i§U 

10 + T J 

* I (18) 


T . 

min 


AT 


min 


Tg sec (p 


2 of 

sf~tAv 




(19) 


for the 183.31 GHz line. 

Table 4 of minimum detectable optical depths is obtained from 
equations 18 and 19. The same limiting values of <*/, yp as before’ 
are used. 

Figures 6 and 7 are plotted with the limiting atmospheric 
optical depth values from Table 2 and the minimum detectable optical 
depths presented in Table 4. The figures show that a radiometer tuned 
to the 22.235 GHz line can measure water vapor up to 85 km under 
optimum conditions with 70 km being a more realistic limit. If, 


Table 4: Minimum Detectable Optical Depths for Passive Resonant Absorption { T g sec 0 ) 
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however, the experiment is performed at 183.31 GHz, the measure- 
ment can be made all the way up to 100 km for the optimum conditions. 
However, this type of experiment should have an observation time of 
10 seconds or shorter due to the relative movement of the source. 

Thus the upper limits of measurement would be less. 


6.2 Active Resonant Absorption 

This is the most versatile method. Global coverage is possible 
although an additional vehicle is required. The source can be a fixed 
frequency broad band source or a narrow band tuneable source. This 
method has the slight disadvantage that permission to opera.te the 
source must be obtained. However, the ability to control the source 
characteristics and the fact that this is a non-thermal source and thus 
will have less noise inherent in its output, make this experiment 
extremely attractive. 

The minimum detectable optical depths for this case are 


T . 

mm 


AP . 
mm 

P c sec tf> 

D 


a?k v 

"Tr” 



for the 22.235 GHz line and 


t . 
mm 


^^min = gk Va v 
PgSOCtf. NTT" 



( 20 ) 


( 21 ) 


for the 183.31 GHz line. 

In the above expressions P is given by equation 15 and has a 

o 

- 14 

value of 8. 84 x 10 Watts. From equations 20 and 21, the following 


table of minimum optical depths is obtained for both lines. 


Table 5: Minimum Detectable Optical Depths for Active Absorption ( sec 
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From tables 2 and 5, figures 8 and 9 are plotted. Figure 8 
shows that if a 22,235 GHz radiometer is used having a signal to noise 
ratio of 10 and an absorption time of 1 second, then water vapor can 
be measured up to 90 km if the true mixing ratio is 10 ppm. However, 
if the observation time is 100 seconds and the signal to noise ratio of 
1 is sufficient, thenwa-ter vapor can be measured above 100 km if the 
mixing ratio is 10 ppm and up to 90 km if the mixing ratio is 0. 1 ppm. 
Figure 9 on the other hand shows that a great improvement results if 
the 183.31 GHz line is used. For a radiometer having a signal '■o noise 
ratio of 10 and an observation time of only 1 second, water vapor can 
be measured up to 85 km if the mixing ratio is . 1 ppm and above 100 km 
if the mixing ratio is 10 ppm. A measurement made with a radiometer 
having a signal ! o noise ratio of 1 and an observation time = 100 seconds, 
can detect and measure water vapor above 100 km even if the mixing ratio 
is only .1 ppm. Thus this measurement seems to be very feasible, 

An experiment having only 1 second for observing time should give 
quite accurate results. The sensitivity of the active experiment can be 
increased by increasing the power of the transmitter. This will then 
move the AP . / P„ curves to the left. However, there is a major 
problem to be considered in this technique. 

At high power levels transitions among the molecular states are 
induced at a rate that is not negligible compared with the collision rate, 
thus invalidating the assumption of thermal equilibrium for the specific 
rotational stite (Karplus et al. , 1948). The absorption line is then 
broadened as the power level is increased with a corresponding decrease 
in the peak resonant cross section. The energy absorbed per unit 
volume will then reach a saturation level and thus perturb the 




Figure 8 22.235 GHz Active Absorption 
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measurement. This means that the level of power used to excite the 

medium may be extremely critical and thus this possibility should be 

examined more closely. This situation is looked at in Appendix D and 

it appears that saturation effects will not be important for stratospheric 

2 

measurements unless the power density exceeds 19 millwatts/m . 


CHAPTER VII 


CONCLUSIONS 

The distribution of water vapor in the upper atmosphere must 
be known if one is to obtain a better understanding of upper atmospheric 
physical processes,, This will be necessary to assess the impact of 
artificially introduced products on the photochemical and radiational 
balances of the upper atmosphere 0 Water vapor may also be respon- 
sible for hydrated ions and noctilucent clouds* Thus a measurement 
of stratospheric and mesospheric water vapor is of considerable 
importance to the scientific community* 

Past measurements of water vapor have been limited to about 
30 km due to the lack of sufficiently sensitive instruments* However, 
due to state of the art advances in low noise receiver development 
(the Schottky barrier diode), weak signals from stratospheric water 
vapor can now be measured* 

The results of this study show that using an occultation 
experiment, water vapor concentrations as low as 0* 1 ppm can be 
detected above 30 km and with sufficient instrument resolution, up to' 
and above 100 km. The 22*235 GHa line of water vapor can only be 
used for an active absorption experiment, if measurements to 90 km 
are desired (see Table 6)» However, the water vapor line at 183*31 GHz 
will enable measurements to be made up to 80 km using the least 
sensitive mode of operation (passive emission). This mode also has the 
advantage that a radiation source is not required and therefore the 
experiment will be greatly simplified and more nearly continuous 
measurements can be made* The active absorption experiment might 
offer the advantage of probing the water vapor line shape more 
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Table 6: Upper Altitude Limits of Water Vapor Measurement 
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accurately, but the emission experiment performed at 183,31 GHz 
appears to be the simplest, least expensive, and most desirable 
experiment. 

The experiment should now be examined using a computer 
simulation. The radiative transfer problem will be examined in much 
more detail especially the effects of combined pressure and Doppler 
broadening on the (signal) line shape. Various atmospheric water 
vapor models will be simulated to determine the relative effect on the 
signal. The orbital and spatial motions of a vehicle will be studied as 
to the restrictions they place on a measurement. The amount of vfrater 
vapor in the viewing path will be varied in order to determine the 
degree of sensitivity of the signal to various regions along the viewing 
path. The measurement should also be modeled to see how far off the 
line center a measurement can still be made. Comparisons between 
several vehicles and various types of measurements will then be 
drawn. This study will provide the theoretical background on which 
to base future experiments. 



APPENDIX A 


WATER VAPOR ABSORPTION COEFFICIENTS 

The absorption coefficient for the pressure broadened 22.235 
GHz line is given by (Croom, 1965} 

2 

K (v, Av, T, N) = 1.05xl0" 28 exp (-644/T). 
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(A. I) 


where 

-3 

N Is the number density of water vapor molecules in cm 

v is the frequency in Hz 

and 

T is the kinetic temperature in °K 

The absorption coefficient for the pressure broadened .183*31 
GHz line is given as (Croom, 1965) 
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(A. 2) 


In previous atmospheric work (Groom, 1965) absorption 
calculations have been simplified by assuming a pressure broadened 
line shape for all altitudes. This has the effect of assuming the same 
line shape for Doppler broadening as for pressure broadening. Above 
70 km where Doppler broadening becomes dominant, this is in error 
since Doppler broadening has a different line shape (Longbothum, 1974). 
The error, however, appears to be slight according to Croom (1965) 
and thus the only account of Doppler broadening in the calculations is 
through the line width constant Av. When both the effects of Doppler 
broadening and pressure broadening are appreciable (above 70 km), 
the value of Av is given approximately by (Croom, 1965). 

A. = (Av p 2 + Av d V /2 (A. 3) 


The pressure broadening is given by (Croom, 1965) 
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(A. 4) 


where 

P is the total atmospheric pressure in mb 

_3 

p is the density of the water vapor in gm m 
T is the kinetic temperature in °K 

The Doppler broadening is given by (Croom, 1965). 
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where M is the molecular weight of the gas. 

For water vapor (M = 18) the expression for Doppler 
broadening is 


Av 


8.45 x 10" 8 Vt* v q 


Hz 


(A. 6) 



APPENDIX B 


DERIVATION OF THE EQUATION OF SENSITIVITY FOR A 
MICROWAVE RADIOMETER VIEWING A THERMAL SOURCE 

Consider the total power superheterodyne receiver (see figure 10) 
monitoring a thermal noise source such as the solar disc. For this 
case A T, the input signal will be broadband random noise. The two 
sided power spectral density of V^ will be (Tiuri, 1964) 

( j u) = | k (T a + T r + AT) = + AT) 

(B. 1) 


The noise input then will be band limited by the bandwidth of 
the tuned RF and IF amplifiers. The bandwidth of the IF amplifier is 
generally smaller and thus has the dominant effect. The bandwidth is 
given by (Burdic, 1968) 


B HF = 


j* |A{jto)l 2 df 

- o 



(B.2) 


where w = 2ir f 


After the RF and IF amplifiers, the spectral density (power) 
is obtained from equation (1) by multiplying it with the power amplifi- 
cation ratio j A(j w) | 


Thus 



Figure 10 Total Power Superheterodyne Receiver 
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<L (jco) = $k( T sn + AT) | A (J'eo) | 2 

(B. 3) 

If we assume that the input signal is' gausslan white noise with zero 
mean and autocorrelation function 


R n-' T » ■ (B. 4) 

Then the output from the square law detector will be given by (Thomas, 
1969) 
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From this 
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The first term is the D. C. r,. "ponent resulting from the 
rectified signal and system noise. The second term is the convolution 
of the input spectrum. 

The D. C. tdrm constitutes the signal power that is measured. 
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The average D. C. power delivered to a one ohm resistor in the 
low pass filter is given by 
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(B. 8) 
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(B. 10) 


At the first stage of the low pass filter, a voltage 
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is injected to cancel the D. C. component due to system noise. The 
difference is a measure of the signal noise temperature. 
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V s E [V ] 
e L e J 


= £k,AT f |A(jw)| 2 d£ 


{B. 12) 


This represents the signal going into the integrator. The low pass 
filter will remove most of the random fluctuations. The variance 
(mean squared normalized power) of the noise fluctuations is given by- 
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The post detect! on bandwidth is given by (Carlson, 1968) 
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Since this bandwidth is much smaller than u may be set equal 

to zero in equation (13). The average D. C. V ^ is given by 
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If equation (13) Is divided by this relation (equation 15), then the 
relative deviation is obtained. 
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For any physical system, there exists an uncertainty principle 
between bandwidth and rtsetime given by (Carlson, 1968) 
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Thus for a perfect integrator with integration time Tj, 
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and 
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(B. 20) 


The sensitivity of the radiometer is defined as the signal which 
will give a D. C. output voltage equal to or greater than the effective 
value of the output fluctuations due to the system noise (Tiuri, 1964). 
Thus if we consider an output signal to noise voltage ratio a, then the 
equation for sensitivity is 


A T 


mm 


"SN 


4 AT 


mm 


a j 1 
V B HF t i 


(B. 21) 


This can be solved for A T, 

MIN 


a T 


AT 


SN" 


Min 


V 


e hf t i - “ 
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There is one problem with the total power receiver. The 
receiver can not distinguish between the change in signal power and 
the change in receiver gain. Thus it is very important to stabilize 
the gain of the receiver. The best receiver for this purpose has been 
that of Dicke (1946) in which he introduced a modulation of signal 
input in order to eliminate the fluctuations induced by receiver 
instabilities. The input of the receiver is switched between the 
antenna and a reference noise power. The rate of switching is at 
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such a frequency that the amplitude of gain fluctuations is negligible. 
The amplitude of the signal then may be determined by means of a 
coherent detector which is driven at the same switching frequency. 
However since the signal is connected to the receiver only half of the 
time, the sensitivity of the Dicke receiver is only one half as high as 
the total power receiver thus 
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For practical bandwidths this will reduce to 
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APPENDIX C 


DERIVATION OF THE EQUATION OF SENSITIVITY FOR A 
MICROWAVE RADIOMETER VIEWING A NON- THERMAL SOURCE 

Consider the total power superheterodyne receiver (see figure 
10) monitoring a continuous wave signal from an oscillator. For this 
case the signal input to the receiver will be 


x(t) = s (t) + n(t) 


(C.l) 


where s(t) = A cos (w o t + 9) and n(t) represents the random back- 
ground and receiver noise. A is a constant and 0 is a random 
variable uniformly distributed over the interval O < 0 < 2 -it. For 
wide sense stationary noise and signal, the autocorrelation functions 
are (Thomas, 1969) 
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From these the power spectral density can be obtained 
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After the RF and IF amplifiers, the spectral densities are 

obtained from equations (4) and (S) by multiplying them with the power 

2 

amplification ratio | A ( j to ) | . 

$ IF = TT k ^SN I w {C. 6) 


and 


^IF Snai (J “ ) = ~ 2 ^ '~ [ 6 ( w " w 0 ) + 6 (w + w Q )] 1 A( jco ) | 2 

(C.7) 

Then the autocorrelation functions at the input of the square law detector 
are 
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If we assume that the input noise is gaussian with zero mean and auto- 
correlation function as given above, then the output from the square 
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law detector will be 
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where the interaction of the noise with itself is (Thomas, 1969) 
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where 
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The interaction of the noise with the signal is given by (Thomas, 


1969) 
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where 
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The interaction of the signal with itself is given by 
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The signal at the input of the detector has been modified by its passage 
through the RF and IF amplifiers. Thus the signal at the input of the 
square law detector is given by 


s(t) = A | A ( j t*>) j cos ^co Q t + 0 + $(to)j 


= A 


j cos [co Q t + v] 


(C. 16) 


where y is still uniformly distributed over the interval o < y <2tt . 
Then the autocorrelation function of the detector output is 
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Thus the square law detector output autocorrelation function is 
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The spectral density at the output of the square law detector 
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Then 
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Thus the square law detector output spectral density becomes 
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The first term in equation (23) is the D. C. component resulting 
from the rectified signal and system noise. The average value of this 
term is 
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If a voltage V o is injected into the first stage of the low pass filter 
where 



(C. 25) 


then the D. C. voltage due to noise temperature will be canceled and 
the difference V is a measure of the signal temperature where 
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The receiver output spectral density is given by 
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However, the low pass filter will pass only those frequency components 
in the vicinity of zero. Thus the output spectral density is 
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or using 
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one can write 
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Another definition of noise bandwidth is (Thomas, 1969) 
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This can be used to write 
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and thus equation (30) can be written as 


<ut n 0“> 


„ ,4 .2 2 

9 f a 2 cr 4cr or ^ 

' ! H <j “) I { 2 " °s 6< “> + * 5 ^: + B° H “ } 


HF 
(C. 33) 


The output autocorrelation function is then 
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The average or D. G» value of the output is 
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The mean square value of the output is 
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The output variance is 
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The relative deviation of the output is given by (Tiuri, 1964) 
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where 
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The sensitivity of a radiometer is defined as the input signal 
which will give a D. C. output voltage equal to or greater than the 
effective value of the outpiit fluctuations due to system noise (Tiuri, 
1964). Thus if we consider an output signal to noise voltage ratio a 
(a figure of merit for the receiver), the equation governing the 
Sensitivity of a receiver is 
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or from equation (38) 
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For a perfect integrator (see Appendix B), the low frequency 
bandwidth can be written 
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Using this relation, equation (41) can be written as 
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This equation has the solution 
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and using equation (32) 
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Thus the sensitivity equation becomes 
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where A now represents the smallest power that can be detected for a 


given output signal to noise ratio. 
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Then since 


AP . 


min 


A 


2 


(C. 48) 


is the smallest power that can be detected at the receiver input, the 
equation of sensitivity for the active source is 
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For a Dicke receiver, the sensitivity will be decreased by a 
factor of 1/2 (see Appendix B), thus 
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In general this will reduce to 


AP 


min 


2t sn ak B HF 


V 


e hf t i 



(C. 51) 



APPENDIX D 


POWER SATURATION OF MOLECULAR STATES 

The expression for an absorption coefficient which is broadened 
by saturation is (Carter and Smjth, 1948) 
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The expression for saturated line half width (at half intensity) is 
(Carter and Smith, 1948) 
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The expression for the dipole moment of the 22,235 GHz 
water vapor transition is (King et al, , 1947) 
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The Einstein coefficient for induced emission then becomes 
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From equation (D.2) saturation effects will be negligible when 


AP << ( A v ) 2 
' o 

Thus the arbitrary criterion will be 
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From the Information in Table 1 the following table of 


maximum power densities is obtained 
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h 

(Av q ) 2 


P 

max 


30 km 

14 

16.0 x 10 

Hz 2 

4.5 x 10 2 

watts 

2 

m 

50 km 

7.3 x 10 12 

Hz 2 

2.0 x 10° 

watts 

2 

m 

70 km 

5. 3 x 10 i0 

Hz 2 

1.5 x 10" 2 

watts 

2 

rn 

90 km 

6.8 x 10 7 8 * 

Hz 2 

1.9 x 10“ 4 

watts 

2 

m 


The cross sectional area illuminated by an antenna of area 

2 

10 m at a distance of 1000 km (the minimum distance between 

7 2 

source and receiver in an occultation experiment) is 1.9 x 10 m . 

Thus the power of the source can be on the order of 3 Kw before 
saturation must be considered at 90 km. At lower altitude* the 

source power can be much greater. 
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